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Abstract 

Recently it has been shown that surface plasmons supported by an interface between a 3+1 
dimensional topological insulator and a metal or between a 3+1 dimensional topological insulator 
with residual bulk charge carriers and vacuum have a polarization that is rotated by an angle 
of order of the fine structure constant a compared to their topological trivial counterparts. In 
this work we generalize this analysis to the more realistic case of thin films, taking into account 
the effect of multiple reflections. In the symmetric case of a thin film surrounded by the same 
material on both sides the polarization of both allowed surface plasmon modes is unchanged 
from the single interface case, even though their dispersion relation is altered by the interactions 
between the two surface layers. In the general asymmetric case the angle is affected as well. We 
give a simple analytic expression for the resonance condition that determines the angle and solve 
it perturbatively for the case of a thin film between two media with almost identical permittivity. 

1 Introduction 

Recently the properties of surface plasmons propagating at the interfaces involving three dimen- 
sional topological insulators (TIs) have been analyzed in [1 , following a similar earlier discussion 
of the two dimensional case in [2J. To allow a non-trivial surface plasmon at least one of the 
two materials has to have a negative permittivity e(w) as would, for example, be found in a free 
Drude metal. One scenario in which this can be realized is an interface between a TI and an 
ordinary metal. A second scenario described in [Tj is an interface between an ordinary insulator 
(e.g. vacuum) and a topological "insulator" with a residual bulk charge carrier density (a sit- 
uation that experimentally is much easier to achieve than a genuine TI). The properties of the 
surface plasmon are very similar in the two cases. 

The analysis in [I used the TPs description in terms of their low energy effective theory. In 
the presence of a time reversal breaking perturbation, such as an external magnetic field, TIs 
are believed to be described by a topological field theory containing an E ■ B term [3]. Such a 
term has interesting effects on the propagation of the surface plasmons. While their dispersion 
relation is only modified at order a 2 (where a denotes the fine structure constant), the surface 
plasmon's polarization is already affected at order a. In particular, while a surface plasmon 
mode on an ordinary metal/insulator interface is purely polarized in the transverse magnetic 
(TM) direction, it was found that for a topologically non-trivial interface the polarization is 
rotated by a non-vanishing angle v v of order a into the transverse electric (TE) direction. This 
effect is very similar to the rotation in polarization experienced by light transmitted through 
or reflected by a Tl/vacuum interface (the Faraday and Kerr angles [3]) and, like the latter, is 
potentially accessible experimentally. 
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The analysis in [1 was performed only for an interface between two semi-infinite regions filled 
by two different materials. Of course in reality TIs will only occupy finite regions of space. In a 
finite geometry such as for example a thin film, the effects of multiple reflections can potentially 
significantly alter the properties of electromagnetic waves interacting with the interface. For 
example, for the case of the Faraday and Kerr reflections it was found in [4] that the thin film 
geometry allows for several new effects. With this motivation in mind, we analyze the properties 
of a surface plasmon propagating on a thin film with two topologically non-trivially interfaces. 
Our main finding is that the rotation of the plasmon's polarization (the angle v p mentioned 
above) remains completely unaffected for a thin film sandwiched in between one and the same 
medium on both sides, while it changes if the upper medium and the lower medium are different. 

The organization of this work is as follows: in section 2 we will review the properties of 
surface plasmons in thin films for the case of standard insulators. In section 3 we will then 
generalize these findings to the case of topologically non-trivially interfaces. 



2 Review of surface plasmons on a thin film 

The basic conditions for having collective excitations of surface modes, the quanta of which 
are known as surface plasmons, has been reviewed in [5]. In particular, we are interested in 
p-polarized surface plasmons (transverse magnetic modes) which are practically accessible. We 
require the plasmons to oscillate tangent to the surface and decay exponentially perpendicular 
to it from both sides: 

Hi = e y H^e l ^ x - ut h lk i z z > 

(1) 

Hn = e y H^e< k " x -^ e~ lk ^ z z < 

and: 

Ei = ^(e x fef - e z kf y(^x^t)^ lZ z > 

o (2) 



Where: fcf = in\ and fcfj = Note that ki and ku need to be real and positive to guarantee 
an exponential decay in the normal direction. These are solutions to Maxwell's equations which 
support surface plasmons provided that the boundary conditions on the components of electric 
and magnetic fields are satisfied at the interface z — 0. This leads to the following condition: 

^ + ^ = (3) 
ei en 

This implies that in order to get surface excitations of this type, one of the media should have 
a negative permittivity. 

The other possible solution to Maxwell's equations is the s-polarized (transverse electric) 
mode. Then one can show that the necessary condition for having surface plasmons is: 

^ + ^=0 (4) 
Mi Mil 

As it is hard to achieve a negative magnetic permittivity in practice, we ignore TE modes from 
now on by assuming real positive magnetic permittivity for both media equal to that of the 
vacuum. 

Now let us consider a three medium configuration with two interfaces at z = and z = —d. 
We call media in z > 0, — d < z < and z < —d regions medium I, II and III respectively. 
The surface plasmon condition is changed due to having two sets of boundary conditions at two 
interfaces. 

One can proceed to satisfy boundary conditions for solutions to the Maxwell's equations as 
given above. However, a simpler way is to look at the reflectivity of an incoming light beam 
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and find out the condition for resonance of the system. At the z = interface the electric field 
for p-polarized incoming, reflected and transmitted light are: 

= ^F(-extf - e z kf)e l ^ x ~^e-^ z z > 
E{ = 4p(e x fcf - e z kf z> o (5) 

En = ^^(-exkfi ~ e^)^^-^ e~ ik ^ z z<0 

Note that for the surface plasmon we are again interested in kf = in\ with real and positive k. 
In this case requiring a solution that has no exponentially growing piece at z > means that 
we are looking for a solution with only reflected, no incoming wave. Even though these names 
are somewhat misleading in this case as the waves aren't propagating in the z-direction, this is 
still mathematically equivalent to asking for a resonance where the total reflection coefficient, 
including the effect of multiple reflections, goes to infinity. 

Continuity of E\\ and if || provide two equations from which we can read the ratio of the 
reflected amplitude as well as transmitted amplitude to the incoming amplitude. To proceed, 
we need only the reflection coefficient from medium II: 

Ei = euk? qfcfi 

Ei e n kf + ejfcf! V ; 

To find the total reflectivity of the system, we should add up the first reflected ray at the 
upper boundary and the subsequent rays which have partially reflected from the lower and upper 
boundaries and then transmitted to the medium I afterwards. For a symmetric configuration, 
when an infinitely long thin layer of thickness d with permittivity en is placed in a medium 
with permittivity ej, the total reflectivity i? 12 of the system can be obtained by summing the 
geometric series corresponding to multiple reflections: 

R = r 12 + r 21 e- 2 ^ d 
12 1 - r 21 r 2l e- 2 ^ d 1 ' 

A resonance would occur when the denominator of total reflectivity meets its zero; this gives 
the dispersion relations of the symmetric layer system [6] : 

eifcfj + enfc I z tanh(^) = 

(8) 

eifcfj + e n fcf coth(^) = 

Either one of these dispersion relations describes an allowed mode. Not surprisingly, we get two 
surface plasmons, one for each interface. In the d — > oo, we restore two copies of the surface 
plasmon condition for two media (jSJ. For small d they mix with each other, which is reflected 
in the modified dispersion relation. 



3 Surface plasmons on a thin film topological insulator 

In the case of a topological insulator while the Maxwell's equations remain unchanged, the 
constituent relations are modified. In SI units : 

D = eE - eoaf(c B) 

(9) 

where cq is the vacuum speed of light and a is the electromagnetic structure constant. The 
parameter 9 can take values (topologically trivial medium) or tt (topologically non-trivial 
medium). 
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As it is shown in [T], the surface plasmons at the interface of a topological insulator and a 
normal medium are not only TM polarized, but they develop a non-zero TE component such 
that at 0(a) the ratio of the TE to TM amplitude is proportional to a ( ej ~ ejl ) . So we need 
to write the most general form of surface plasmons' electric and magnetic fields as a linear 
combination of TM and TE modes: 



Ei = -e y E$ + ^(e x kf - e z kf)e l ^ x ~^e lk ^ z z > 
E n = -e y E u + ^^(-e^fi - e z k u )e i{ - k ^ x - ut ">e- ik n z z<0 
Using the relation: = — £ x (cB) gives: 

Bi = e y ^ + ^-{e x kf - e z kf)e^ x - ut h ik " z z>0 

Bu = e a |n + ^-(-e x k u - ^kf^e 1 ^-^ e~ lk " z z < 



(10) 



(11) 



Using the constituent relations ([§]) in both media, we can easily build H and D fields. Then 
boundary conditions on the field components at z — give the dispersion relation as well as TE 
to TM amplitude ratio. The former as is shown in [I] is only modified at 0(a 2 ) while the later 
gets modification at 0(a): 

. . Ett a AO cuen „, ,aA6 , . 

tan t/ p =-^- = ~ - ~ 7T " \ + ° — ( 12 ) 
En 7r c (ei-e n ) it 

where: Ad = 9\ — 0jj and we should keep in mind that e\e\\ < and: ei + en < PQ. 

We aim to see how the plasmon polarization angle v p will change in case of a thin layer of 
topological insulator placed in a normal medium. Similar to the non-topological thin film, we 
analyze the problem by looking at the resonance state of a light beam reflecting from the film. 
The incoming, reflected and transmitted electric and magnetic fields at the boundary z = are 
respectively: 

Ei = -e y E { l + ^(-e x kf - e z kf ) e i(fc?*-"t) e -ik?* z>0 
E[ = -e y E^' + ^{e x k z -e z kf)e t( - k i x ~^e lk i z z>0 (13) 
E u = ~e y E u + E ^(-e x kf l - e^y^n*""*^-*" 2 z<0 

and: 

B{ = e y ^ + ^-(-e x kf - e z kf)e i( - k i x -^e- lk ~i z z>0 
Bi = e y ^ + ^-{e x kl~e z kf)e t{k ^ x - ul ^e lk ^ z z>0 (14) 

Bu = e y ^ + ^-(-e^fcfj - e z kf l )e t ^ x ~ ul ^e- lk " z z < 

Using the constituent relations ©, we can also form H and D fields. Boundary conditions 
are the continuity of Eu and Hu as well as continuity of D± and B± at z = 0. This gives six 
equations, out of which two are redundant. For any of the continuity equations to hold for all x 
and t we need k x to be equal to kfj. So we are left with four independent linear equations with 
six unknowns: E°, E h Ef ', E{, E u and E u . 

Treating E® and E\, that is the amplitudes of the two polarization components of the incom- 
ing wave, as input parameters, we can solve for the corresponding amplitudes of the reflected 
light in terms of a reflection matrix R defined as: 

f)=m). as) 
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That is the reflected amplitudes E®' and E[ are linear combinations of the incoming ones, E® 
and E\\ if one sends a purely TM (or TE) polarized light toward a topological insulator, the 
reflected light has mixed TM and TE components. Explicitly we find 



0/ _ -gfHiiCfcff+^fc^+fc^ 

1 ~~ eil(fcf) 2 +ei(fcf I ) 2 +/cffef I (£i+eii+£ (^) 2 ) 

E , _ gi [^(^)fc?fcn1+gih(fcn 2 -^i(fcn) 2 +^fcn(-^i+^n+^(^) 2 )] 
1 ~~ £ll(fe I 2 ) 2 +ci(fc I z I ) 2 +'c I z fe I J! I (ei+eil+«o(^) 2 ) 



(16) 



where we have set /ii = fin = fio- 

For the symmetric case of a thin TI film immersed in a medium (so that e\ = em and 
similar for 9 and /x), it is most convenient to first diagonalize the matrix R. An eigenvector of R 
corresponds to a polarization for which the incoming light gets reflected off the interface without 
being rotated. In terms of these eigenvectors we can then perform the resonance analysis as in 
the topological trivial thin film. We'll find two different resonance conditions for the two different 
eigenvectors. One of them is the rotated TM mode we are interested in. That is, its polarization 
only differs from TM by order a and it needs negative e in one of the media to be supported 
by the interface. The second polarization is the rotated TE mode. Its resonance condition can 
not be met unless one of the /x's is negative and so we'll disregard it as before. Naming the 
eigenvalues of the matrix i? 12 as (r"o) 12 and (r), 2 , the well-defined reflection coefficients then 
are: 



Oo)i2 



-fcffclW^(^) 4 +2e ( ei + e „)(^)3 + (e 



£ll(fc I i ) 2 +ci(fef I ) 2 +fc-fej ! I ( e i+e I i+eo( 2 ^ £ ) 2 ) 



, , , = • -^l(fcn) 2 +fcffc I ' i yi 2 ('^)4+ 2eQ ( ei +e II )(^) 2 + ( e i- 

In order to impose two resonance conditions: 



(17) 



l-(ro) 2 i(ro) 2 ie- 2K « d = 
1 - (rhi(?hie- 2Kud = 



(18) 



we need (r~o) 21 and (r) 21 coefficients instead. These are eigenvalues of matrix i?2i obtained 
from (|16|) after interchanging I and II indices. These conditions give two dispersion relations. 
Defining the variable x = ri- , it is not hard to see that these relations are a polynomial equation 

quartic in x and has been modified by terms of at least 0(a 2 ) similar to the result obtained 
in [1] for the topologically non-trivial two medium system. Now let us assume that medium 
I is a normal insulator such as vacuum, and medium II is a topological insulator, then en is 
necessarily negative and also: |en| > ei as stated in Ref [3]. Then one can show that the only 
solutions for the variable x which reduce to ((5J) at A8 — come from the second equation in 
(|18p . These are the rotated TM modes we are looking for; the other resonance condition would 
correspond to the rotated TE modes. 

In order to extract the angle v p we need to look at the ratio of the TE to the TM amplitude 
for the rotated TM mode. To do so, we should look at the corresponding eigenvector of the 
matrix R21 as defined above and apply the dispersion relation to its components. But one can 
easily check that these eigenvectors are independent of x, and we have: 



tonfo) s — = ,.„.„, „A8. (M) 



CO 

„2\ 



Expanding equation (fT§|) up to 0{a 2 ) gives: 

to W = "f;& + 0((f) 2 ) *><fc (20) 
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which shows that tan(f p ) will not change in case of a thin film topological insulator placed in 
a normal medium compared to a semi-infinite topological insulator in proximity to that normal 
medium. 

This result seems to be a consequence of the symmetric configuration we have considered. 
If the upper and lower layers are filled with different substances, we should expect a different 
polarization angle. To prove this claim, we proceed as follows. The general resonance condition 
should be written in the matrix form: 



12x2 - ^21-^236 2l%d — 



(21) 



Note that R21 and R23 can, in general, not be diagonalized simultaneously, so we need to 
diagonalize the matrix M = 12x2 — R2iR-23e~ 2Kd as a whole and finding the dispersion relations 
by setting its eigenvalues equal to zero. Although this can be done numerically quite easily, a 
concise analytic solution for the general asymmetric case can not be obtained. But in order to 
show how the polarization angle is affected in this configuration analytically, we adopt a simple 
asymmetric case where the top layer is vacuum and the bottom layer is only slightly different 
from vacuum while the middle layer is a topological insulator: 



em 



ei = e = 1 

e n = e 
e (1 + 5) = 1 



(22) 



With 6 <C 1. And we set cq = 1 and c\\ = c. Also: 



dr = x + S Xl 



(23) 



where xq — yo is the solution to the symmetric resonance equation. We can plug these values 
into the eigenvalues and eigenvectors of matrix M and expand them at order 5. Although the 
eigenvalues, which determine the dispersion relation after being set equal to zero, depend on 
x\ and 2/1 explicitly, the eigenvectors on the other hand are independent of x\ and y\ at this 
order and we only need to use the solution to xq = yo of the symmetric case, which we found 
before. Note that there are two acceptable solutions for Xq which reduce to the two modes of 
© at A9 = 0. That is they correspond to the expected two rotated TM modes supported by 
the two interfaces. One can show that there are two corresponding polarization angle at O (aS) 
as following: 



tan(i/jp) = — —-fj 



(l-e) 



tan(^) = -^^ 



1 + (l-«) - { l+e^r+e(-l+e^f d + U K 6 ) 
1 _ (l-e) -(- 1+ e^f+c(l+e^f d + U K 6 ) 



0(a 2 ) 
0(a 2 ) 



(24) 



where e < 1 and |en| > ei- 

This clearly shows that if we place a thin film topological insulator between two different 
media, the polarization angle v v changes compared with a configuration where a thin film 
topological insulator is placed in one medium. 
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